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1  INTRODUCTION 

The  random  walk  is  a  classical  mathematical  problem.  The  definitive  statement  of 
the  problem  was  given,  many  years  ago,  by  K.  Pearson' .  Pearson  required  to  know  the 
probability  of  finding  a  random  walker  (or  -  to  paraphrase  Ref  2  -  "a  drunken  man  in 
open  country")  within  a  given  radius  of  his  starting  point  after  N  equal  steps.  The 
random  walker  changed  his  direction  in  a  random  manner  after  each  step. 

In  the  classical  random  walk  all  angles  through  which  the  walker  may  turn  after 

each  step  are  equally  probable.  The  angle  probability  density  function  is  uniform.  In 

this  paper  the  non-uniform  (or  non-isotropic)  walk  is  considered  in  which  the  step  angle 

probability  density  function  is  not  uniform.  Probability  density  functions  which  are 
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mapped  on  to  the  unit  circle  have  been  discussed  by,  for  example,  Wintner  ,  who  gives 
further  references.  In  this  paper  the  existence  of  phase  probability  density  functions 
is  taken  for  granted;  they  can  be  easily  justified  on  physical  grounds.  Consider  the 
example  of  an  electromagnetic  wave  incident  upon  a  point  in  space.  The  wave  is  scat¬ 
tered.  Suppose  now  that  the  phase  of  the  scattered  wave  is  measured  relative  to  the 
unscattered  wave  at  the  point  of  emission.  Suppose,  further,  that  random  fluctuations 
in  the  path  length  between  the  scatterer  and  the  measurement  point  make  the  phase 
measurement  uncertain,  then  the  phase  has  a  probability  density  function.  If  the 
fluctuations  are  such  that  the  path  length  has  a  Gaussian  density  then  the  resulting 
phase  probability  density  function  will  be  a  Gaussian  with  its  'tails'  periodically 
folded  into  the  interval  -  tt  to  +n  ,  because  measurement  of  phase  is  only  unambiguous 
within  this  interval.  This  case  has  been  examined  by  Wintner  in  Ref  3  (it  corresponds 
to  the  case  of  the  parameter  1  in  equation  (1)  of  Ref  3  equal  to  2).  It  is  examined 
further  in  section  4  of  this  paper  but  from  a  physical  viewpoint. 

R.  von  Mises  seems  to  have  been  among  the  first*  to  consider  the  idea  of  a 
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probability  density  function  mapped  on  to  a  circle  and  von  Mises  density  function 
(Ref  4,  equation  (12)  and  Fig  5,  page  494)  is  employed  in  this  work  in  section  5.3. 

In  connection  with  Ref  4,  see  the  comments  by  Greenwood  and  Durand  in  Ref  6,  section  2. 
Greenwood  and  Durand  call  von  Mises  density  function  the  'circular  normal  distribution' 
but  it  would  be  more  apt  to  call  the  Gaussian  (or  'normal')  density  mapped  onto  a  unit 
circle  the  circular  normal  density,  although  in  this  paper  that  function  will  be  called 
the  theta  density  for  reasons  which  will  become  apparent  in  section  4. 

The  non-isotropic  two-dimens iona 1  random  walk  appears  to  have  been  considered 
first  by  Rayleigh^  who  considered  the  walk  in  which  the  angle  density  functions  of  each 
step  were  uniform  over  an  interval  not  equal  to  2nr. 

Returning  now  to  the  example  of  the  f luctuating-point  scatterer,  consider  an 
experiment  in  which  the  phase  and  amplitude  of  the  scattered  wave  is  measured  at  a 
number  of  intervals  of  time.  It  is  assumed  that  sufficient  time  is  allowed  to  lapse 
between  each  measurement  so  that  the  samples  are  not  significantly  correlated.  The 
measured  vectors  are  added  together  on  the  complex  nlane  to  form  a  random  walk.  It  is  the 


*  Although  see  footnote  (1)  on  page  8  of  Ref  5  where  it  appears  that  P.  L£vy  had  con¬ 
sidered  probability  distributions  mapped  on  to  the  unit  circle  several  years  earlier. 
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object  of  the  work  presented  in  this  paper  to  investigate  the  phase  probability  density 
function  of  the  resultant  of  such  a  random  walk.  In  particular,  the  case  of  most 
interest  was  that  in  which  the  phase  density  function  of  each  step  was  almost  uniform 
and  it  is  to  this  case  that  most  attention  is  directed.  It  will  be  shown  that,  for  a 
large  number  of  steps,  even  slight  deviations  from  uniformity  in  each  step  phase  density 
function  can  make  thu  phase  density  of  the  resultant  sharply  peaked.  In  other  words  the 
effect  of  the  random  walk  is  to  amplify  the  non-uniformity  in  the  step  phase  probability 
density  function. 

One  avenue  of  approach  to  the  problem,  at  least  for  a  large  number  of  steps,  is  to 
invoke  the  central  limit  theorem  (see  Ref  8,  §5-4).  The  central  limit  theorem  predicts 
that  for  a  large  number  of  steps  the  real  and  imaginary  components  of  the  resultant  of 
a  walk  on  the  complex  plane  will  each  tend  to  have  a  Gaussian  distribution.  One  would 
then  work  out  the  first  and  second  moments,  write  down  the  joint  density  function 
(Ref  8,  §8-2),  and  then  compute  the  phase  density,  using  the  standard  methods  of  trans¬ 
forming  probability  functions  given  in,  for  example,  Ref  8,  §3-6.  This  method  was  not 
adopted  for  two  main  reasons.  Firstly  it  was  desired  to  see  how  far  one  could  go  in 
constructing  an  analysis  valid  for  a  small  number  of  steps  as  well  as  a  large  number. 
Secondly  an  analysis  based  on  the  central  limit  theorem  rapidly  becomes  unwieldy  and 
awkward  to  handle  in  the  case  of  the  phase  density  function.  The  analysis  which  follows 
is  based  largely  on  the  theory  of  generalised  functions  (concentrated  on  smooth  hyper¬ 
dimensional  manifolds)  given  by  Gel'fand  and  Shilov  in  Ref  9.  No  attempt  has  been  made 
to  be  consistently  rigorous  and  some  of  the  analysis  is  based  on  physical  ideas.  In  the 
final  section  the  results  of  this  analysis  are  checked  against  numerically  computed  phase 
probability  density  functions  obtained  by  a  Monte-Carlo  method. 

Throughout  this  paper  the  term  'probability  density  function'  is  used  in  preference 
to  'probability  distribution'  in  order  to  avoid  confusion  with  the  completely  different 
usage  of  the  term  distribution  in  generalised  function  theory. 


2  THE  JOINT  PROBABILITY  DENSITY  FUNCTION 


2.1  Introduction 


In  this  section,  the  joint  probability  density  function  (which  henceforth  will  be 
called  'the  density  function')  of  magnitude  and  phase  of  the  resultant  of  a  two- 
dimensional  non-uniform  walk  is  derived.  It  is  first  introduced  in  an  intuitive  way, 
and  then  in  sections  2.2  to  2.4  it  is  obtained  in  a  more  rigorous  manner. 


The  density  function  is  denoted  by  p(a,6)  and 
resultant  of  a  random  walk  within  a  distance  da  and 

£q 

ae  on  the  complex  plane  is: 


dp  «  pia.Cldad0 

Since  it  is  a  certainty  that  the  vector  has  somewhere 


the  probability  of  finding  the 
within  an  angle  d>*  of  the  point 

n 

•  -  .c.uicu _ _ _  _____ 

- (>) - - 

r>- . 

on  the  complex  plane:  l  t;t  i  on/  _ 
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p( a, 6) dad 8  =  1 


(2) 


2.1.1  A  one  step  'walk' 

Consider  a  'walk'  consisting  of  one  step  of  length  r  in  a  direction  <t>  ,  <t 
being  measured  anticlockwise  from  the  real  line  to  the  vector  jr  ,  then  the  ’probability 
density  function'  corresponding  to  this  'walk'  of  fixed,  certain,  length  r  and  fixed, 
certain,  angle  <f  is 


p(a,6)  =  6(a  -  r)6(B  -  <f)  (3) 

p(a,9)  is  to  be  regarded  as  a  generalised  function  defined  in  terms  of  the  delta  func¬ 
tions  6(a  -  r)  and  6(9  —  <p) .  It  is  necessary  that  r  >  0  .  The  angle  delta  function 
is  periodic,  with  a  period  of  2tt  .  Thus 


p  (a,  0)d9da 


0  “IT 


CD 

'll 

0  -IT 


6(a  -  r)6(9  -  (J>)drd<)>  =  1 


(4) 


for  all  e  e  and  r  >  0 

Suppose  now  that  $  has  a  continuous  density  function  over  -it, it  .  Let  this 
function  be  P(<1>).  Then 


p(a,9)  =  P(P) 6 (a  -  r)  -  5(a  -  r)  j  P(*)6(6  -  4>)d<j>  , 

-IT 


(5) 


where  a  well-known  convolution  property  of  the  delta  function  has  been  used,  see  Ref  9, 
§5.2,  equation  (3),  page  104,  So  (5)  is  obtained  by  convolving  (3)  with  the  angle 
density  function.  Note  that  this  is  easily  extended  to  the  case  where  r  is  also  con¬ 
tinuously  distributed  for  then: 


p(a,0) 


00 

/ 


Q(r)6(a  -  r)dr 


/ 

-TT 


P(t)6(0  - 


for  r  >  0,  9  G  [~  nt7r'] 


(6) 


2.1.2  The  two  step  walk 

Now  suppose  that  the  walk  has  two  steps  of  fixed  and  certain  length  Tj  and  r^ 
at  fixed  certain  angles  of  d>  1  and  p?  .  Then  evidently 

/  i<f>l  i$2  \ 

p(a,9)  «  5(a  -  l^e  +  r2e  j  (9  -  tp)  (7) 

m 
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where  the  resultant  of  vectors  r^  and  r_ ^  is 


resultant 


I".6 


+  r26 


i'll 


(8) 


If  <J>  j  and  <f>  are  allowed  to  have  any  value  with  probability  densities  Pj(f-j)  and 
P2 (^2^  then  convolving  (7)  with  P  j  ( ( )  and  P^ves: 


p(a,9) 


7T  TT 

'  // 


P.^l 


)P2(«2)f 


(*' 


i<i, 


V 


r2e 


5(9  - 


iti)d5jdi); 


(9) 


and  again  the  result  may  be  extended  to  the  case  where  r^  and  have  continuous 

densities  Q((rj)  and  Q?(r by  convolving  over  Tj  and  .  Physically,  (9)  may 
be  regarded  as  the  sum  of  all  combinations  of  vectors  of  lengths  r(  and  r^  over 
all  angles  weighted  by  their  respective  probabilities,  the  integrations  being  regarded 
as  limits  of  sums  over  the  probability  densities.  It  will  be  seen  further  that  if 
<t>  j  and  d> 2  are  not  independent  then  Pi  ^  j  )?2^2^  *s  rePlaced  by  the  j°int  density 
p(<t>j<}>0)  in  equation  (9),  although  in  this  paper  only  independent  steps  are  considered. 

2.1.3  The  N  step  walk 

It  is  straightforward  to  continue  the  argument  to  the  random  walk  of  N  steps 
each  step  having  a  length  r^  ,  k  G  [l,N]  and  angle  <S  ,  G  [-it, it]  .  Then 


1 

where 


N  . 

p(a,9)  =  5^a  -  |  ^  r^e  k!ji(°  - 

k=1 


'll)  , 


the  resultant 


_  l 


IN 

z 


1*.  •  . 
r  e  k  !elU 
k 


(10) 
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If  the  vectors  have  continuous  phase  angle  density  functions  C d>^ ) ,  then 


p(a,9) 


TT  Ti 

hi 


VV 


W5 


(“ 


5(8  -  0)d$ 


d(fv 


02) 


where  there  are  N  convolutions  over  the  and  the  <t>  are  all  independent.  Also 

K  K 

rk  >  0  and  $k  G  [-ir.ir]. 


(12)  is  the  required  expression  for  the  joint  (a, 9)  density  function  of  an  N 
step  non-uniform  walk.  In  general  two  of  the  integrations  may  be  performed  (since  there 
are  two  delta  functions)  but  the  other  integrations  are  very  difficult.  The  density 
function  has  been  obtained  here  in  an  intuitive  way  but  it  can  be  shown  to  follow  from 
the  characteristic  function  of  the  walk  and  this  is  the  subject  of  sections  2.2  to  2.4. 
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2.2  The  characteristic  function 

The  angle  density  function  of  each  step  is  periodic  over  -ir,ir  and  is  defined  by 
a  Fourier  series: 


W 


B  (k>e 
n 


03) 


where  the  notation  B^(k)  means  the  nth  Fourier  coefficient  of  the  kth  step,  B  (k)  may 

TT 

be  complex  and  BQ(k)  =  I  since  /  =  1  ' 


The  characteristic  function  «  f(£,s)  and  is  defined  as 


f(C,C)  =  (exp  i(£x  +  ?y>)  , 

see  Ref  8,  section  4-3,  page  52,  also  Ref  19,  section  3.2,  page  41 


(14) 


where 


and 


U 

x(*,  ...  ♦N;r]  ...  rN)  =  ^  rk  cos  *k 

k=  1 

N 

y(*i  ^N’r i  =  £  rk sin  \  • 

k=  I 


(15) 


(16) 


The  average,  implied  in  the  notation  (  )  ,  is  taken  over  the  angles  <t>^  . 

So 


r  r 

<exp  i p  ^  rk  cos 

+  q  rk sin  *y  * 

i 

07) 

N 

(  exp  i(Srk  cos  <{>k 

1 

+  ?rk  sin  4>k)  ) 

(18) 

N 

(exp  i(£rk  cos  4>k  ■* 

■  Crk  sin  *fc) >  . 

(19) 

1 


The  density  function 
over  the  <f>k  gives: 

of 

the  angles  <f>k  is  given  by  (13). 

Hence  performing  the  averaging 

N 

t"t  r 

IT 

00 

«  "1 

fU.c)  -  III 

1 

2w 

J  exp  i(£rk  cos  4>k  +  Crk  sin  <J>k) 

y  B  (k)  exp  in<J>kd<t>k|  .  (20) 

k-1  L 

-IT 

n=-°° 

Put  £  =  a  cos  n  and  C  3  a  sin  n  ,  then  integrating  the  Fourier  series  term  by  term 
gives 

N  00  7T 

g(a,n)  =  j""j'  ^  Wk)  I  exp[iark  cos  +  ni (4>fc  +  n)]  d<t>k  (21) 

k“l  n“-*°  -x 


N  °°  IT 

=  j  |"  ~  ^  B^(k)  exp  n*n  I  exP^“ric  cos  \  + 


(22) 


k=l  n=- 


Now 


ii 

/ 


exp(iar  cos  $,  +  r»i<f>.  ) d4>  »  2xi  Jn(ar  ) 


(23) 


where  J  is  a  Bessel  function  of  order  n.  See  Ref  10,  §2.2  (5). 
n 


Hence 


g( 


o,n)  =  TTE  B^(k)  exp  in(n  +  tt / 2)  Jn(ark)  . 


(24) 


k=  1  n=-°° 


The  sum  over  n  in  (24)  is  absolutely  convergent.  This  may  easily  be  seen  by  putting 

B  (k)  =  1  .  The  sum  is  then 
n 


exp  in(n  +  it/ 2)  Jn(ark)  =  exp(ia  cos  n) 


(25) 


which  is  absolutely  convergent  (Watson  Ref  10,  §2.1).  Since  all  the  B^(k)  ate  finite, 
■ie  less  than  some  constant  K  the  series 


oo 

E*. 


(k)  exp  in(n  +  tt/2)  Jn(ark) 


is  also  absolutely  convergent. 

So  the  series  in  (24)  may  be  multiplied  together  and  the  finite  product  may  be  evaluated. 
Multiplying  out  the  product  in  (24)  gives 


N  <* 


TTE*. 

k=»  1  nc-“ 


(k)  exp  in(n  +  ir/2)  Jn(ark) 


oo  oo 

E  -  Iv‘>-b« 


(N)  exp  im  (n  +  it/ 2)x 


V 


X  Jn  (ar  )  . . .  Jn^ar^  ,  (26) 


where  m  =  n  ^  +  ...  +  n^  . 
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g(a,n)  =  ^  •••  ^  B^O)  ...  B^N)  exp  im(n  +  tt/2  )  Jn^  (ur  t  )  . . .  Jn^ai^)  (27) 

n,=-“  v-" 

which  is  the  required  expression  for  the  characteristic  function  expressed  in  radial  and 
angular  frequencies. 

2.3  The  density  function  as  a  sum  of  Hankel  transforms 

The  density  function  is  given  by  the  inverse  Fourier  transform  of  the  character¬ 
istic  function  (Ref  8,  §4-3).  In  polar  frequency  coordinates  this  is 

go  7T 

p(x,y)  =  — J  J  exp  -  ia(x  cos  n  +  y  sin  n)g(a,n)  dn  ada  (28) 

^  0  -7T 

or  if  x  =  a  cos  6  and  y  =  a  sin  0 

oo  n 

p(a,6)  =  —^2  J  J  exp£-  iaa  cos(n  -  0 )]g(o,  h)  dn  adct  (29) 


Now,  consider 


I  =  j  exP[~  iaa  c°s(n  -  e)]g(a,n)  dr) 


IT  »  00 

■h  j  exP[-  iaa  cos<P  -  9)]  X!  Z!  Bn 


(1)  ...  B^fN)  exp  im(n  +  x/2)  x 


n  _  =—oo  n  ——oo 

1  Ti 


x  Jn^  (ar  j  )  . . .  Jn^iar^)  dn  (31) 


OO  00 

r  -e  Bn(1)  •••  VN)Jnl(url)  ...J^Carj.)  exp  imx/2  x 
V“  V” 

TT 

I  exp[- iaa  cos(n  -  9)  +  iron]  dn 

—  TT 

77 

~  j  exp£-  iaa  cos(n  “  6)  +  imnj  dn  =  exp  im(9  -  rr/2)  Jm(aa) 


CO  oo 

I -I  Bn(l)  ...  Bn(N)Jnj  (at|)  ...  Jn^ar^)  exp  im0  Jm(oa) 


n  ss— co  n  =— « 

1  N 
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and 


p(a,6) 


co 


0 


(35) 


Hence 

p(a,8) 


oo  oo  oo 


—  I  jT  . . .  y  Bn(  1 )  . . .  Bn(N)  Jnt  (ar  J )  . . .  Jn^(arN>  exp  imP  Jm(aa)a  da  ,  (36) 


0  nl=_"  V” 


which  is  a  sum  of  Hankel  transforms.  This  is  the  required  result.  Notice  that  if  the 
angle  density  functions  are  uniform  then  the  only  non-zero  term  in  the  sum  is  that  for 
which  =  ...  =  n^  =  0  .  In  which  case, 


ie 


p(a,6) 


p(a) 


~  I  JQ(ar  ( )  . . .  JQ(arN)  JQ(aa)  a  da  , 

0 

OO 

a  I  Jp(arj)  ...  J0(arN)Jp(aa) a  da  , 

0 


(37) 


(38) 


which  is  Klyver's  integral  for  the  uniform  walk  -  Ref  11  and  Ref  12,  also  Ref  19, 
section  4.4,  page  93. 


Integrals  similar  to  those  in  (36)  have  been  studied  by  several  writers.  The  most 
recent  appears  to  be  Bailey  ^  who  gives  references  to  earlier  work.  Watson'^  also  gives 
further  references.  In  particular  Bailey  shows  that  (Ref  13,  equation  (10.3)): 


OO 

f  Jp(bt)Ju(a  t)Jv(a  t)  ...  J  (a  t)t‘  U  V’"k+'dt  =  0  ,  (39) 

J  I  2  k  n 

0 

when  b  >  a,  +  a„+  . . .  +  a  and  all  a,  are  positive. 

12  n  k 

In  our  case  o  =  u  +v+  ...  +k  and  this  obviously  reduces  to  the  integrals  in  (36). 

The  result  in  (39)  means  that  the  integral  has  bounded  support  -  or  a  finite  span  - 
physically  the  integral  has  a  finite  span  because  the  resultant  a  cannot  be  greater 
than  Tj  +  rj  +  ...  +  rN  •  The  probability  that  a  >  tj  +  ...  +  r^  is  zero  and  this  is 
the  physical  interpretation  of  Bailey's  result.  Unfortunately  Bailey's  technique  cannot 
be  used  to  evaluate  the  integral  when  a  <  Tj  +  ...  +  r^  .  Under  certain  circumstances 
the  resultant  a  may  be  confined  to  an  annulus.  Suppose  that  one  step  r^,  say  is 
greater  than  the  sum  of  all  the  others,  ie  >  r7  +  ...  +  r^  .  Then  the  probability 
that  a  <  r j  -  -  . . .  -  r^  is  zero.  This  is  also  contained  in  Bailey's  result  since 

a  <  Tj  ~  r2  ~  ■’*  ~  rN  **  rj  >  r2  +  r3  +  •••  +  rfj  +  a  anc*  this  is  just  a  matter  of 
redefining  the  walk.  It  will  next  be  shown  how  (36)  may  be  transformed  into  (12). 
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2.4  Transformation  of  the  density  function 
We  have 


p(a,0) 


.  B  (n)Jnj  (ctr  |  )  .  ..Jn^(ar^)  exp  im0  Jm(aa)  a  da  , 

.  (36) 


and  the  Fourier  coefficients  B  (k)  are  given  by 

n 


i  -irL$ 

Bn(k)  =  /  P(4>k)e  d<fk  ‘ 


(40) 


Replacing  the  B^(k)  in  (36)  by  (40)  gives: 


p(a,6) 


oo  77  oo  co 

£  /  fH\)d\  Y,  Z  Jn«<ar»>  JnN(“rN)  * 

0  -tt  n  j  =-“  t\i=_0° 

x  exp[im0  -  i(nj$j  +  ...  +  n^<j>j,)] Jm(aa)  ada  ,  (41) 


where 


IT 

f  P<Vd*k 

—  TT 


means 


N 


repeated  integrals  on  the  same  interval. 


ie 


TT  7T  3t 

j  P(*k)d*k  =  |  ...  J  p«1)...P(+N)d*|...d*lj  . 

-TT  -TT  -TT 


(42) 


i  0 

Now  the  point  on  the  complex  plane  at  which  the  density  function  is  measured  is  £  =  ae 
Subtracting  the  last  step  from  £  gives: 

c  =  a  -  r  (43) 

—  —  -N 


so  that 

See  Fig  1 . 


2  2 

a  +  r„  -  2ar„  cos(P 
N  N 


(44) 


i 

Let  m.=  30  m  ^  and  m^  =  n^  +  .  Then  the  sum  over  in  (41)  may  be 

3  k=l  u1  '* 

written 


Y  Jn^(ar^) Jm(aa)  exp  i(m0  -  n^^) 


co 

exp  itnN_,0  l  +  ®N_))(aa)JnN(arN)  exp  in^e  -  *N> 


V 


cn 

sf 


(45) 
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The  sum  on  the  RHS  in  (45)  may  be  summed  by  a  Bessel  function  addition  theorem  (Graf's 
formula),  see  Watson  Ref  10,  §11.3  (2).  Then 

0» 

^  J(nN  +  mN-!)(°ia)JnN(arN)  6XP  inN(?  "  V  =  JmN-l(ac)  exp  imN-l(n  "  f  )  »  (46) 


where  c  and  n  are  defined  in  Fig  I. 


^  Jn^iar^)  Jm(aa)  exp  i  (m6  -  n^f^)  =  exp  im^^o 


And  so 


p(a,0) 


°°  11  1X1 

j  ada  j  P(bk)d<fk  ^  ...  ^  Jnl(otri) 


0  -V 


exP[imN-|n  “  +  •••  +  1Vt4N-l)3JlVl(ac) 


Evidently  this  process  can  be  continued  by  subtracting  the  next  step  r,,_j  from  £  and 
so  on.  When  the  walk  has  been  reduced  by  N-!  steps  then: 


co  TT  co 

p(a,e)  =  f  ada  jf  P(4>x)d$k  ^  Jn j  (ar  ( )  Jn (  (av)e 


in(  (>-:,) 


0  -it 


n^=-.oo 


1  =  !«  “  IN  -  •••  _  I2!e 


Applying  Graf's  formula  for  the  last  time  gives 


co  TT 

p(a,0)  =  ~  [  ada  j 


P(i )di.  Jn(aw) 
k  k  0 


To  proceed  further  it  is  necessary  to  reverse  the  order  of  integrations  in  (51);  diffi¬ 
culties  then  arise,  however,  because  the  integral  over  a  diverges  at  the  upper  limit, 
In  fact  as  a  -*  10 

aJQ(aw)  ~  /2a I trw  cos  (aw  -  tt / 4 )  .  (! 


It  will  be  shown  that  (51)  can  be  interpreted  as  a  generalised  function, and  to  this  end 
a  function  qt(a,0)  is  introduced: 


f  *7  I 
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CO  ^71 

<’t(a'e)  “  -t  I  3  exp  ("7")  d0  f  p(Va* 


kVaw) 


and  Ura  q  (a,0)  »  p(a,6)  ,  t  is  an  integer.  Hence  p(a,0)  is  the  limit  function  of  a 

t-x»  ^ 

sequence  of  function  qt(a,0)  as  t  -*  .  The  convergence  is  now  dominated  by  the  expon- 

ential  factor  and  the  order  of  integration  is  reversed  to  give* 


Tf  CD 

qt(a,6)  =  £  f  P<Vd*k  /  “  eXp('7f) 


jQ(aw)  da  , 


where  w  ■  |a  -  £  r.  i  . 

~  I 

Now  consider  a  fixed  walk  consisting  of  N  steps  each  of  length  r ,  . .  .  at  fixed 
angles  of  $  . , .  .  Let  the  resultant  of  this  fixed  walk  be  s  exp  ill  . 

s  exp  ill  =  r;  exp  i|,  +  ...  +  r^  exp  i*N 


w  =  a  -  s  , 


x  *  2 

qt(a,0)  »  ~  j,  Pfifi^ld*^  j  a  exp  JQ  (a  I  a  exp  i6  -  s  exp  i\l|)da 


Graf's  formula  is  now  applied  but  in  reverse.  See  Fig 


■^0 (a I—  “  —  “  /  dn(aa) Jn(as)  exp  in(C-  -  i) 


TT  oo  2  ^ 

qt(a,0)  *  j f  P(4i^)d'}i^  J  a  exp  ^-•^p^da  Jn(aa) Jn(cts)  exp  in(fi  -  if) 

~r  o'"  n=-°° 

The  Neumann  series  in  (60)  is  uniformly  convergent,  so  integrating  term  by  term: 

■rr  oo  oo  ^ 

qt(a,6)  “  j  P(<fk)d(j>k  exp  in(6  -  i if)  j  Jn(aa) Jn(as)  exp  ^--^—^ada 


°°  2 

j  Jn(aa)  Jn(as)  exp  ada  =  2tP  exp[-t‘(a“  +  s‘)] In(2t^as) 


where  I  is  the  modified  Bessel  function  of  order  n  ,  see  Watson  Ref  10,  §13.31  (1) 
n 

Also,  for  large  t 


•>  1  2 
I  (2t“as)  exp(2t  as) 

n  2t/?aF 


See  Watson,  Ref  10,  §7.23. 


Jn(aa) Jn(us)  exp 


(-i) 


t  1  2.  .2 

adct  ~  —  — —  exp  -  t  (a  -  s) 

V  7T  /as 


and  this  is  a  delta  convergent  sequence.  See  Gel'fand  and  Shilov,  Ref  9,  §2.5,  ex  2. 

The  integral  J  Jn(aa) Jn(as)ada  can  therefore  be  interpreted  as  a  delta  function  and 

0 

we  have 


Jn(aa)  Jn(as)  ada  =  5(a  -  s)  , 

r  as 


and  this  is  independent  of  n  . 


j  P(4k)dik  a  t(a  -  s>  )  exP  in<s  "  •  <6 


The  Fourier  series  —  \  exp  in(0  -  il)  defines  a  'row  of  deltas'  at  ?  -  +  2nn  , 

n=-°° 

or  in  the  case  of  this  problem  the  delta  function  at  d  -  0  since  the  domain  of  9  -  ip 
the  interval  -ir,ir  . 


p(a,0)  =  i  -  s)5(6  -  <!■)  , 


p(a,{S)  = 


/.../ 


P,  (d; , )  ...  PN(dN)d4>  j  ...  d*N6(a  -  s)6(e  -  <i)  (6 


I r j  exp  id,  +  ...  +  rN  exp  if^l 


s_  =  j  s !  exp  i4' 


which  is  the  required  result. 
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It  has  thus  been  shown  how  (68),  ie  (12)  which  was  obtained  in  an  intuitive  way  in 
section  2.1.3  may  be  obtained  from  the  characteristic  function  of  the  walk. 


Two  results  which  were  shown  above  are  repeated,  that  is 


where 

and 


Jg(ctz)a  da 


-y=  6  (x  -  y)6(a  -  6) 
/xy 


2 

z 


2  2 

x  +  y  -  2xy  cos(a  -  8) 


00 

l  Jn(ax)  Jn(ay)  ada 
0 


y=  «(x  -  y) 
/xy 


(71) 


(72) 


(73) 


In  passing  it  is  pointed  out  that  by  making  use  of  the  convolution  property  of  the  £ 

9 

function  (Gel-fand  and  Shilov  §5.2)  it  is  straightforward  to  show  that  (71)  leads  to 
Neumann's  integral  theorem  (Watson^  §14.3  (1))  and  (73)  leads  to  Hankel's  integral 
transform  theorem  (Watson'^  §14.3  (3)  and  the  rest  of  Chapter  14). 

3  SOLUTION  OF  THE  JOINT  DENSITY  FUNCTION 

In  this  section  an  attempt  is  made  to  solve  (68)  directly.  It  will  become  apparent 
that  by  and  large  (68)  is  intractable  but  it  may  be  integrated  in  one  or  two  cases. 

3. 1  The  two  step  walk 

The  steps  have  lengths  r.j  and  r^  and  angle  density  functions  P ^  (d> ^ )  and 

P2<*2>' 

Then;  „  _ 

p(a,e)  =  J  I  Pj  (4>1)P7(£2)6^a  -  j  r  (e  1  +  r2e  2|^6(0  -  ilJdifjdOj  .  (74) 

“TT  -7T 


Now  let 


I  l<t>I  X<t>2 1  i0 

c  =  ! r  j  e  +  r2e  |e 


See  Fig  3. 
Then 


2  2 

rj  +  r2  -  2r j  r2  cos  8 


2  2 

rf  +  r2  +  2r;r2  cos(<{>2 


since  S  =  rr  —  (<t>  —  <t>  )  ,  see  Fig  3. 

Also 

2  2  2  2  2 
Tj  *  c  +  r2  ~  2cr2  cos  a  “  c  +  r2  ”  2cr2  cos^2  “  ^ 

since  a  «  <j>2  -  iJj 


(75) 


(76) 


(77) 


>6 


-  cos 


[  2  2  2t 


cdc /  r,r2  ,  .  cdc 

/  '  \  l  2V*  J  “<V20 


where  A(r(r9c)  is  the  area  of  triangle  r^c  • 

So  transforming  the  integral  over  in  (74)  to  an  integral  over  c  gives: 


p(a,e) 


l  P2(*2U*2  I  Pi[-2  -  COS~1(  — 2-r--r2  ~ )] 

_Tr  lrl"r?l 

2  r  /r2  -  c 

*  6[9  “  ^2  +  cos~'  (~  2ci 


5 (a  -  b)  x 


2  2  2 


-1  /r|  ~  C  ~  r2  \  cdc 

\  2cr2  /J2i(rir2c) 


(80) 


11 

-JJL—  f 

(r.r2a)  J 


P2(<j>2)d<t2P1  $2  -  cos 


Note  that  in  the  integration  over  $2  there  are  two  delta  functions  in  the  interval 


♦ 2  G  [-*•*] 


And  so 


p(a,8)  = 


which  is  the  required  result. 

Note  that  when  the  angle  density  function  is  uniform  P^  =  P2  *  l/2ir  and  then: 


p“  ‘ 


ir6(r]  r2a) 


which  is  the  result  for  the  uniform  walk  given  by  Rayleigh  and  others. 

3.2  The  three  step  walk 

The  steps  have  lengths  r^,  r2  and  r^  and  angle  density  functions  P  ^  (d>  ^ )  , 
P2(4>2)  and  P^H^)  •  Th®  walk  is  shown  in  Fig  4.  It  will  be  noted  that  the  walk  has 
been  divided  up  into  two  two-step  walks.  If  aj  is  the  resultant  of  rt  and  r2  a 
new  walk  may  be  defined  which  is  the  sum  of  a^  and  r^  .  The  density  function  of  a^ 
is  given  by  (82) : 


r*l 
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p(a, ,?, ) 


4tvp'[’ 


0  +  cos 


X  P2i  I  +  C°S 


W)] 


-l 


In  forming  the  density  function  p(a,6)  it  is  now  necessary  to  convolve  over  a 
a(  is  now  continuously  distributed  over  the  interval  £  [lri  -  r2 I  ’  ri  +  r2 


the  density  function  of  the  walk  formed  by  a^  and  is: 


rl+r2  ttt  /  l  2  l\  ill 

I,  l  !*&?’■  ^ 

r  1  — r2 1  L  ’ 


2  2  2' 


2  2 


p (a, 0 )  = 


r  /  2  2 

„  .  -If  rl  “  al  "  r2  \ 

L'  \  2air2  /J 


x  p2  |e ,  +  cos 


p3(*3)  x 


x  a 


i.  6  i  ^ 

-  la^e  !  +  r3e  ^|^5(6  -  '(OdOjdbjdaj 


The  same  transformation  as  before  is  now  applied: 


Let 


so 


also 


if*  i  &  , 

l  I  3 |  ll 

=  'a,e  +  r3e  !e 


c‘  =  aj  +  r3  +  2a j r3  cosifj  -  S j )  , 


2  2  7 

at  =  c  +  r3  -  2r 3c  cos(<t>3  -  1) 


and 


Also 


1  =  $3  -  cos 


d9  j  =  2cdc/A(a j r3c) 


M  ■ 


-i 


6,  =  1-  -  cos 


Hence 


(84) 

since 
.  Hence 


(85) 

(86) 

(87) 

(88) 

(89) 

(90) 

(91) 
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r  1  +r2  IT  a!+r3  „  ,  /  Z  2  Zv  /  2  2  2  . 

(^M.  -  -"(^)  •  -"P&O]  - 

X  P3($3)6(a  -  c ) 6 j^O  -  *3  ♦  cos-(a_l)]^ 


2  2  2, 


2  2  2  , 


p(a,0)  = 


-  cos 


- \  d40  da 

,r3c)  3  1 


Performing  the  integrations  over  c  and  4*3  gives: 

ri+C2  ,  r  /  2  2  2  \ 

p  (a,  9)  =  f  ^  ^  p,  Lcos-' 

]  A(r  r  a  )  Afa.^a)  \  2ar3  ) 

iVr2 ' 

-i/'!'V'iy  -./rr*?-'2\  -i/arrrr2\] 

’“s  1-TvT-JJ 


x  P^  0  +  cos 


(93) 


And  this  is  as  far  as  it  is  possible  to  go  unless  the  angle  density  functions  are  uniform 
in  which  case  the  integral  can  be  performed  and  leads  to  a  complete  elliptic  integral  - 
see  for  example  Ref  14  and  also  Ref  15  Appendix  b.3. 


3 . 3  The  N  step  walk 


(2  9  9  v 

rl  ”  a  r9  \1 
- ~ -  I  I  which  is 

2at2  )i 


In  (82)  the  angle  density  functions  are  P^l  $  +  cos 


1  /  2  2  2k  /  2  _  2  _  2»  - 

P  +  cos  — 2lF2 -  J  ~  cos  ( - 27jTy — y 

and  this  is  P^fO  +  (4,,  -  0)  -  (42  ”  4  j )]  or  P  ^  C  4  j )  .  Likewise  in  (93)  it  will  be  seen 
that  the  angle  density  functions  are  P|[9  +  (<$>3  -  0)  -  ( ct  3  -  4  )  +  ( 4  ^  -  4 ,)]  or 
P ! (4 ! )  and  so  on. 

Hence  it  will  be  seen  that  extending  the  analysis  to  four  steps  from  three  is  a  matter  of 
adding  another  triangle  to  the  walk  and  another  angle  density  function,  ic 


ai+r3  ri+r2 


P(a,n ) 


I  I  Mr,^)  A(a|r3a2)  Afa^a)  P I I  )P2(V  V4  3)P4(C4)  da  !  da2 


|arr3l  ir,-r2 


(94) 
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where 

*  1 

i 

<> 

■F" 

-  *3> 

*  ( v  2  ”  $3)  “  ( &  2  ~  ^  ^  * 

(95) 

$  =  0  + 
v  2 

($,  -  $)  -  ($, 

4  4 

-  *3> 

+  ^2  ”  ^  * 

(96) 

»  + 

”  3 

1 

w 

1 

O 

-  S) 

9 

(97) 

+ 

ii 

<r 

(*4  "  ^  ’ 

(98) 

where  the  angle  differences  are  written  in  terms  of  the  lengths  of  each  triangle  in  the  walk. 

It  is  thus  possible  to  extend  the  result  to  N  steps,  although  the  integrals  appear 
to  be  totally  intractable.  It  is  not  possible  to  find  the  phase  density  function  by 
evaluating  the  joint  density  function  p(a,0)  and  then  integrating  over  a  .  The  prob¬ 
lem  of  finding  p(G)  is  pursued  further  in  section  5  of  this  Report,  but  in  the  next 
section  the  most  interesting  of  the  angle  densities  -  the  Gaussian  density  -  is  examined. 

4  THE  GAUSSIAN  PHASE  DENSITY  FUNCTION  OR  THE  THETA  DENSITY 


Consider  the  following  physical  situation:  a  wave  is  incident  upon  a  point  in 

space  and  is  scattered.  The  phase  of  the  scattered  wave  is  measured  at  some  point.  The 

position  of  the  point  is  uncertain.  Suppose  that  the  distance  between  the  point  and  the 
phase  measuring  observer  has  a  probability  density  function  which  is  Gaussian.  The 

measurement  of  the  phase  has  an  ambiguity  of  2nr  ,  n  being  an  integer.  Hence  phase 

angles  greater  than  t  are  mapped  back  on  to  the  interval  — *  to  n  as  are  phase  angles 
less  than  -t  .  The  phase  angle  density  function  as  measured  by  the  observer  therefore 
consists  of  a  Gaussian  density  with  its  'tails'  periodically  folded  back  onto  • 

The  function  may  be  represented  on  the  interval  by  a  Gaussian  periodically 

repeated  at  intervals  of  2r  .  This  is  illustrated  in  Fig  5.  The  Gaussian  phase  densi  ty 
function  is  then: 


P(0) 


e  e  .  (99) 


Evidently 


/ 


p(o)de  = 


I  i  4 


de 


o>  2ti 


/ 


exp 


(100) 


The  function  P(0)  will  be  recognised  as  a  theta  function  with  an  imaginary  period.  It 
is  in  fact  the  third  of  Jacobi's  four  theta  functions  and  in  the  notation  of  Whittaker  and 
Watson  (Ref  16,  Chapter  21) 

-T 
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P(0)  =  t»  (6/2,  io2/2ir)/o/2¥ 


(101) 


The  corresponding  Fourier  series  may  be  written  down  immediately  (Poisson's  summation 
formula,  Ref  16,  Chapter  6,  Example  18,  page  124,  also  section  21.51,  page  474) 


P(e) 


i  £  exp(-  4-) exp  nie 

n=-°° 


(102) 


P(0)  will  henceforth  be  referred  to  as  the  theta  density,  this  term  being  considerably 
shorter  than  Vlaussian  phase  density'.  It  has  been  assumed  that  the  theta  density  is 
'centred'  on  8=0.  Centring  the  density  function  on  6^  merely  adds  a  phase  shift 
on  to  each  Fourier  coefficient  of  course: 


P(6) 


n=-« 


exp  ni8 


(103) 


The  theta  function  is  very  rapidly  convergent.  Rewriting  (102): 

2. 

p(e>  =  27  t1  +  2  exp(~  ' 

let 


2 

P(0)  =  —■  [1  +  2  exp  cos  9  +  2  exp(-  2c2^  cos  26  + 


(104) 


a  /  2  =  it 


Then 


P(9)  =  —  [  1  +  0.08643  cos  8  +  6.9 


_6 

10  cos 


26  + 


(105) 


Hence  to  a  very  good  approximation  when  o  /2  >  tt  ,  P(6)  can  be  approximated  by  the  mean 
term  and  first  harmonic  only.  This  is  the  case  which  is  of  interest,  as  mentioned  in  the 
Introduction  and  this  approximation  is  used  later. 

In  passing  it  is  mentioned  that  P(6)  may  be  written  in  a  number  of  different 
forms.  If  the  cosines  of  multiple  angles  in  (104)  are  written  as  polynomials  in  cos  0 
it  is  possible  to  prove  the  following: 


P(6)  *  ■—  cos  0  +|  cos2 6  *  -i  (t?|“  +  tfj)  cos30  +~  (i?^V  +  4t3^3)  cos**  6  +  .  .  .  } 


.  (106) 

where  Roman  numerals  denote  derivatives  with  respect  to  z  and  the  first  and  fourth  theta 
functior 
example 


functions  ^  and  and  derivatives  are  evaluated  at  z  =  0,  q  =  e  2a  ,  ie  for 


*1  =  3Z  VZ,q)| 


z=0,  q-e 


-2c 


(107) 


The  proof  is  straightforward  but  somewhat  tedius  and  will  not  be  given  here. 
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5  PHASE  DENSITY  FUNCTION 

Returning  to  equation  (12): 


IT 

p(a,6)  =  ^ 

TT  / 

..  [  ?,(♦,)  ...  pN<y«(a-  1 

N  i<j>  \ 

V* rke  Ms(8  -  il)d<t>!  . 

.  .  ddi 

-*TT 

\ 

Integrating  over  a  gives: 

■n  7i 

l  P,  (♦j)  ...  Pn(«n)6(6  -  d*j  ...  cHn  (108) 

-IT  -n 


and  this  is  the  required  formula  for  p(8).  The  integrations,  however,  are  not  at  all 
straightforward.  Firstly,  attention  is  directed  towards  the  delta  function  6(9  -  ip) . 

5. !  5(3  -  ») 

The  angle  is  the  phase  of  the  resultant  of  a  fixed  walk  with  step  lengths  r^ 
and  step  phases  of  $.  •  That  is: 

&  =  arg  (r  |  exp  i<t  +  . . .  +  exp  i$N)  .  (109) 

If  ij»  is  limited  to  a  half  plane  (any  half  plane),  then 

arg  (r(  exp  i$  +  ...  +  rN  exp  i$N)  =  tan 

Outside  the  half  plane  there  is  ambiguity  because: 


r  sin  if .  +  ...  +  r„  sin  t}>  - 

— - -p - - - —  .  (110) 

r,  cos  $,  +  ...+  rN  cos  4>N 

m 


r  j  sin  4>  j  +  • 

•  +  rN  Sin  V 

-1 

=  tan 

'r (  sin(ij> j 

+  rnt)  +  . 

•  +  rN  sin(*N 

+  nt)  " 

Tj  cos  |  +  . 

•  +  rN  cos  y 

Tj  CO  S  C  j 

+  nn)  +  . 

.  +  rN  cosO, 

+  ntr) 

J 

for  any  positive  or  negative  integer  n  . 


(Ill) 


Hence  if  SO  -  v)  is  written  as 


5(9  -  *) 


sin  <J>j  + 
cos  <1 1  + 


”•  +  rN  Sin  *n\] 
+  rN  COS  *n/J 


(112) 


we  have  the  required  delta  function  at  s|j  =  9  and,  in  addition,  'spurious'  delta  func¬ 
tions  at  i(i  =  9  ±  nn  .  There  are  then  two  delta  functions  in  the  interval  •  The 

effects  of  the  additional  delta  are  considered  in  section  5.1.2. 


5.1.1  Transformation  of  5(9  -  Q 


We  have 


r !  sin  S  j  + 
tj  cos  I |  + 


sin  i p 
cos  5 


N_ 

N 


tan 


(113) 
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Now  change  the  variables,  so  that  Uj  ■  <f>j  -  8  etc.  This  is  a  simple  linear  rota¬ 
tion  by  an  angle  8  with  no  change  of  scale  and  does  not  scale  the  6  function.  Note 
that  in  (108)  the  convolution  is  performed  by  rotating  the  delta  function  relative  to  the 
plane  containing  the  various  combinations  of  vectors  which  make  up  the  walk.  The  trans¬ 
formation  simply  means  that  the  delta  function  is  now  held  stationary  and  the  planes  con¬ 
taining  the  combinations  of  vectors  are  rotated  relative  to  it. 


tan  i);  = 


r  sin(4>  s  +  8)  + 
Tj  cos(<J>j  +  9)  + 


tan(ijj  -  9)  = 


r i  sm  u |  + 
r^  cos  Uj  + 


+  rN  sin(4>N  +  6) 
+  rN  COs(<t>N  +  6) 


+  rN  Sln  UN 
+  rN  COS  “N 


And  so 


6(0  -  if.)  =  6 


[  \r]  cos  u,  ♦  ...  ♦  rN  cos  ^  ) j 


Generalised  functions  are  uniquely  determined  by  their  local  properties.  Indeed  in 
showing  that  a  delta  function  5[f(u<...  u^)]  is  equal  to  another  generalised  function 
k(u,  ...  u^dfgiuj  ...  uN)],  sa'’,  it  is  only  necessary  to  consider  the  behaviour  of 
f(Uj  . . .  u^)  and  g(u(  .  .  .  i.  a  sufficiently  small  neighbourhood  of  the  support  of 

the  delta  function.  The  support  of  the  delta  function  is  the  set  £a^  for  which 
f(uj  ...  i^)  =  0  ,  or  g(u(  .  uN>  =  0  .  In  other  words  the  support  is  the  N  -  I  dimen¬ 
sional  hypersurface  (m  nifoiu)  f(Uj  .  . .  i^)  =  0  .  In  connection  with  those  local  proper- 

t  .9 

ties  of  generalised  functions  which  are  assumed  here,  see  Gel  fand  and  Shilov  ,  §A1.3, 

_j  /  r.  sin  u  ■*  . . .  +  rN  sin  \ 

pages  144-146.  Note  that  in  (116)  tan  (  - : -  }  is  smooth  and 

\  r,  coSU)  +  ...  +rNsinuN; 

infinitely  differentiable  in  a  'sufficiently  small  neighbourhood'  of  the  manifold 


Tj  sin  Uj  + 


+  r^  sin  Ujj  =  0  .  Hence  expanding  tan 


sin  Uj  + 

■j  COS  Uj  + 


+  rN  sin  Un' 
+  rN  COS  V 


about  the  manifold  r,  sin  u,  + 


+  rN  Sln  UN 


0  by  Taylor's  series  gives: 


tan  [f(u 


where 


tan"  [f (a,  ...  a^] 


t\ — 1 — -vi  (u  - 

...  y]23U*  Ju=a  k 


+  0(u  -  a)* 


n 

3f  (s 

‘  L  ^ 


f(U]  ...  Ujj)  = 


aN^(uk  -  a^)  +  0 (u  -  a)' 


r.  sit  u.  +  ...  +  r„  sin  u„ 
i  1  N  N 


Tj  COS  U^  + 


+  rN  C0S  UN 
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Expanding  f(u^  ...  u^)  about  the  same  manifold  obviously  gives: 


N 

f(ul  "  V  =  Z  (3'  •"  aN)(uk  "  V  +  °(u  “  3)2 

k=l  k 


Hence  in  the  neighbourhood  of  the  manifold 

-I 


tan 


[f(ul  ...  Ujj)]  =  f(Uj  ...  uN)  +  0(u  -  a) “  . 


( 1  20) 


(12)) 


Hence  everywhere, 

fiftan” 1 { f (u j  ...  uN)|]  =  5[f(U]  ...  u„)]  .  (122) 

Now,  if  Q  is  a  non-vanishing  function  of  Uj  ...  and  P  is  a  function  of 
Uj  .  .  .  u^j  which  can  be  zero,  then: 

6(PQ)  =  Q~'<5(P)  (123) 

9 

see  Gel'fand  and  Shilov  ,  §17,  page  236.  This  can  also  be  easily  proved  by  a  Taylor 

series  expansion  about  the  manifold  P  =  0.  Note  that  r,  sin  u,  +  ...  +  r„  sin  u%,  can 

1  1  N  N 

be  zero,  but  (r j  cos  u(  +  ...  +  r^  cos  u^)  can  never  be  zero,  in  fact: 


! r ,  COS  u,  +  ...  +  rN  cos  u^!  <  r,  +  ...  +  rN 


(124) 


Hence 

<[ 


rr,  sin  u,  +  . . .  +  r„  sm  u,,n  .  . 

— - !— - — ^ =  (r,  cosu.  +  ...+  r  cos  u  )«(r  sin  u  +  ...+  r  „sin^) 

r ,  cos  u,  +  . . .  +  r„  cos  u„T  11  N  n  i  ' 

LI  1  N  NJ 


(125) 


And  so: 


5Ln-/rjsinu4.-:  +  r.N-8inM] 
L  \r|  cos  u)  +  •••  +  rNCOSUN/J 


(r  |  cos  u  j  +  . . .  +  cos  u^)  6 (r  j  sin  Uj  +  ...+  rNsinu^,) 

.  (126) 

and  in  this  form  the  delta  function  is  amenable  to  calculation. 

5.1.2  Some  properties  of  5(9  -  if) 

There  are  two  deltas  For  each  set  satisfying  r^  sinu^  +  ...+  r^  sin  u^,  =  0  . 

One  delta  lies  on  the  positive  half  of  the  real  line  and  one  lies  on  the  negative  half. 

This  can  be  rectified  by  replacing  (r  ^  cos  u^  +  . .  .  +  r^  cos  u^)S(r  sin  u^  +  . . .  +  r^  sin  u^,) 
by: 


[“*(£  rk cxp  n,0]  "  i[ 


rk cot  v  *  L  'k  cos 

k- 1 


N 

\  6  £  rk  sin  \  (,27) 

J  Lk=  i 


and  then  on  the  positive  half  of  the  real  time 


)  r,  cos  u,  ■  =  V  r,  cos  u,  ,  and  on  the 
/  k  k  ,  k  k 

k  k 


24 


I l  rk  cos  uk! 

k 


l  rk  C°S  "k 

k 


So  the  deltas  on  the  negative  half  of 


negative  half 

the  real  line  subtract  leaving  only  those  on  the  positive  half.  Formula  (127)  therefore 
enables  the  analysis  to  be  extended  over  the  whole  plane,  and  although  the  modulus  term 
makes  the  expression  awkward  to  use,  its  presence  is  inevitable. 

So  equation  (108)  can  be  written: 


p(9) 


TT 

TT 

0 

r  A 

N 

/- 

.  /  P^Uj  +  0)  .. 

••  V“x  * «  1 

*  )  rk  cos  V  + 

)  rk  COS  Uk 

-7T 

-7T 

k=l 

k^1  -1 

r  n  , 

[Z  rx  Sin  UkJdU'  "■  dUN 


028) 


where  u 


4>l  -  6  etc,  and  all  the  P^iu^  +  are  periodic  with  period  2ir 


Now  replace  the  +  9)  by  their  Fourier  series  representations: 


VUk  +  e)  -  2 


V  Y  vk)  exp  inkfuk +  o) 

v~°° 


029) 


Hence: 


N  ”  00  IT  IT 

p(e)  =  ^  Yj  ’•  Z  Bn0)  •••  Bn(N)  exp  ime  [■■■  I  £Xp  1  (Vl  +  •••  +  VV 

nl=-“  nN=-“  -71 

x  i(l  £  rk cos  +  Y  cos  Uk)5(Y  Ck  sin  “k) 


rk  s*n  \  du  j  ...  du^,  .  (  1  30) 


where,  as  usual,  m  =  ^  +  ...  +  n^  .  It  is  evident  that  if  m  is  odd  adding  *  on  to 
the  Uj  •  •  •  ujj  reverses  the  sign  of  exp  i(njU(  +  ...  +  n^u^).  Now  the  delta  on  the  nega¬ 
tive  half  of  the  real  line  generated  by  the  term  cos  d^r^sin  u^  is  negative 

k  k 

and  that  on  the  positive  half  of  the  line  is  positive.  Therefore,  if  m  is  odd 

£  r^  cos  u^5(j|r^sin  picks  out  all  combinations  twice  since  both  deltas  cover  the 

k  k 

same  interval,  the  changes  in  sign  cancelling  each  other.  If  m  is  even  then  adding  a 
if  on  to  the  U|  .  .  .  u^  leaves  the  sign  of  exp  Kn^u^  +  ...  +  n^u^,)  unchanged.  In  this 

case  the  deltas  generated  by  )  r^  cos  >^*(1  r  sin  Uk)  cancel,  and  the  result  is  zero. 

k  k 

Therefore  £  r,  cos  u, &()  r  )sin  u.  pick  out  the  odd  harmonics  in  p(P).  A  similar 

k  k  k  'k  k/  k 

argument  for  |  £  r^  cos  u^  |  d  (][  sin  u^)  can  be  made,  but  with  a  sign  reversed, 
k  k 


And  so 


generates  the  odd  numbered  harmonics  in  p(G)  ,  (131) 


Sr‘ “■  siouk) 

(i  £rk cos  v)6(|\ sin  \) 


generates  the  even  numbered  harmonics  in  p(G) 


Unfortunately  the  modulus  in  (132)  renders  its  application  somewhat  awkward.  There  is, 
however,  one  class  of  problems  for  which  (131)  on  its  own  is  sufficient  and  these  are 
problems  involving  phase  density  functions  which  are  non-zero  only  on  an  acute  segment. 

5 . 2  The  phase  density  function  on  an  acute  segment 

If  the  P^(^^)  are  limited  to  an  acute  segment  and  are  zero  outside  this  segment, 
then  the  resultant  of  the  walk  must  also  lie  within  this  segment.  This  is  also  true  if 
the  segment  is  an  open  half  plane.  It  is  also  true  if  the  segment  on  which  the  P^(‘  ) 
are  non-zero  is  a  closed  half  plane,  with  the  exception  of  one  open  half  plane,  for  then 

the  resultant  still  lies  in  the  open  half  plane. 

Suppose  that  P  (}^)  =  0  for  >  a  ,  a  <  -’/2  ,  then  p(0)  =  0  for  if-1  >  a  . 

Now  if  p ( P )  is  computed  on  the  basis  of  formula  (131),  ie  if  p ( 6 )  is  written: 


=  (£)  Z  Z  B"(l)  •••  Bn(N)  exp  ime 


n  ,  n  ”’cu 


/■■•/ 

-IT  -7T 


exp  i(njUj  +  . .  .  + 


^  c°s  ^(Z- sin 


rk  cos  \N  Z,rk  Sln  \ 


Then  only  odd  harmonics  are  obtained,  and  a  function  of  the  type  shown  in  Fig  6  is 
obtained.  For  -a  <  0  <  a  the  ’true'  p(S)  is  obtained  and  for  -r<e<a-r  and 
it  -  a  <  3  <  ir  ,  -  p(P)  is  obtained.  Of  course,  since  it  is  necessary  that  p(B)  >  0  onl? 
the  result  for  -a  <  P  <  a  has  any  significance,  but  withn.  this  range  it  is  correct, 
and  p(0)  must  be  identically  zero  outside  it. 


On  the  other  hand  if  (  1  ]  r.  cos  u  !  )5  (  /  r,  sin  u,  \  had  been  used  in  (133)  then  the 

V  k  k  k  /  \  k  k  V 

function  shown  in  Fig  7  would  have  been  obtained,  :’c  p(0)  for  -  a  <  C-  <  a  and  p(-  ) 

again  for  -  n  <  0  <  a-r  and  r  -  a  <  0  <  it  . 

It  is  straightforward  to  show  that  a  function  such  as  is  shown  in  Fig  6  where  for 
-  a  <  E<  <  a  the  function  equals  p(G)  and  for  nit-a^3^n1,  +  a,  n  odd,  -nf^)  has  a 
Fourier  series  which  must  have  only  odd  harmonics,  ic  if 


r,  sin 
k 


exp  inO 
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Then 


=.  ■  ['*<->"*']/ 


I  +  (- 1  )  ‘  ‘  |  |  f (6)  exp  -  in®  d®  , 
-a 


and  for  a  function  of  the  type  shown  in  Fig  7  there  are  only  even  harmonics: 


+  (-  l)nJ  I  f (6)  exp  -  inS  d6 


One  of  the  many  sequences  which  define  a  delta  function  is 


S(x)  = 


1  im  I 


/ 


exp  i  ix  d;  , 


see  Gel ' fand  and  Shilov  ,  §2.5,  page  38  (5), 
And  so: 


^ (r |  sin  U|  +  . . .  +  rN  sin  uN>  =  ^  j;  f  exp  { "(r|  sin  u 


rN  sin  UN>  d 


Applying  this  expression  to  (133)  gives: 


CC  CO 


o(9) 


=  (^)  Z-ZB-o)  ••■vn)  exp  iw" 


n,=-  n.j=-,c 


I  •••  f  exP  i(V,  +  •••  +  vV [ 


r,  cos  u.  +  ...  +  r  cos  u  1 
N  A  J 


1  ^  “1 


-TT  -TT 


*  lim  —  j  exp  i^(r^  sin  +  ...  +  r^,  sin  u^)  du^  ...  du^. 


or,  rearranging: 


p  (0  ) 


(1)  ...  B  (N)  exp  im11 
n 


n  •— 00 
N 


N  7T 

II -I 

k=i  -n  -tt 


r,  cos  u,  exp(ir,u,  +  ir  ,  sin  u,) 
k  k  It  1  1 


exp(iryiN  +  ir^,-  sin  u^,) 


du  ...  du^. 


(135) 


(135) 


(137) 


■  .  (138) 


(I  30) 


.  (140) 
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Now  consider  the  integral: 


I  =  I  cos  u.  exp(in,  u  +  ir,  •'  sin  u,  )  du, 
j  k  k  k  k  K 


There  are  two  cases  j  ^  k  and  j  =  k.  Both  are  covered  by: 


I  =  j"  cos  01  exp(ir,  •"  sin  u  +  in^u)  du 


2-(-)n+"‘i[j  ^  ( r r  )  +  J  (;r  )] 
>-  n+ »  k  n-a  k  J 


Hence  for  .1  =  0  and  a  =  1  : 


I  explin^u 


A 


ir^-;  sin  u)  du  =  (-)  2-'Jr^(rv %)  , 


cos  u 


nk+1  \ 

exp(in^u  +  i  sin  u)  du  =  (-)  2r-jp-r  Jn^  (r^ ') 

k“ 


where  a  standard  Bessel  function  recurrence  formula  has  been  used  (Watson'*'',  §2. 
page  17  ( 1 ) ) . 

Therefore : 


j  ■■■  f  rk  cos  uR  exp(in)u)  +  ir  ( c  sin  u()  •  •  •  expfir^.u^  +  ir^f  sin  t^) 


du , 


=  (2-)S(-)W+>  -p-  d n j  ( r  j i )  . .  .  Jn^.  (r^f) 


where,  as  usual,  m  =  n|  +  ...  +  n^. 
Hence 


,1 


p(  ‘)  *  —  1  im  j  d 


I  - 1  ■. 


( 1 )  ...  D  (N)  exp  imB 
n 


nr~"  v 


(-) 


m+  1 


7  '’"AA  •••  A(rN^  Y, 


\ 


(142) 


(143) 


(144) 


1  2 , 


A’ 


(145) 


(146) 
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p(0) 


-  £  t  •••  £  ».«>-•.» 


(N)  (-)m+  'm  exp  imMim  [  Jn  ^  (r  (  r  )  .  .  .  Jn,.  ( rv r )  4- 


n.=-“ 


V 


(147) 


If  m  is  even  then  Jn  (r  ■')  .  .  .  Jn^,(rN,r )  is  an  even  function  and  Jn^r^")  .  .  .  Jn^  (r,. ■.)  / 
is  an  odd  function.  The  integral  is  then  zero.  Hence  there  are  only  odd  numbered  har¬ 
monics  in  (147). 

Equation  (147)  can  alternatively  be  written: 


P(“)  =  lim  /  ~ 


■V  L  k~  1  1  n_ 


(-)  Rn(k)  exn  in^Un^tr,  r) 


I]  ■ 


(148) 


and  this  is  a  useful  form. 

;  \ 

The  series  >  (-)  B  (k)  exp  irt  vJn^(r^t)  is  a  Neumann  series  and  the  coefficients 

V“* 

B^(k)  are  the  coefficients  in  the  Fourier  expansion  of  the  angle  densi ty  func tion  of  the 

kth  step.  For  example,  if  each  step  is  certain  to  be  on  the  real  line  at  r  =  0  then 

the  density  is  a  delta  function  and  then  B  (k)  =  1  . 

n 

Hence 


QC  CC 

Y  (_)\<k)  exp  irvjnk(v:)  =  Y (_>rv  exp  v-w 

nk=_o°  v* 


r)  =  exp  -  ir^'  sin 


(149) 


see  Watson'^*,  §2.22,  page  22. 
And  so 


p(6)  =  -rr  lim  [  i 


y2  i  jr-  exp  -  iisin  r^  =  ~  H™  J 


s  cos  0  exp  -  isf  sin  2  dl 


.  (150) 


where  s  =  )  r, 


p(0)  =  s  cos  05(s  sin  0)  =  cos  96(sin  ?) 

=  l  (“)n6(9  “  mr) 


(151) 


ignoring  all  delta  functions  except  that  within  the  valid  domain  of  9  ,  i<:  that  at 
3=0,  leaves: 

p(0)  =  6(0)  ,  (152) 
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which  is  the  correct  result,  for  the  resultant  is  certain  to  lie  along  the  real  axis  at 
■•=(>. 

This  is  a  particularly  simple  example  and  unfortunately  there  do  not  appear  to  be  any 
more  step  angle  probability  density  functions  which  are  zero  outside  an  acute  segment, 
which  lead  to  a  series  which  can  be  summed  to  give  a  simple  function.  For  example  if 
the  step  density  is  uniform  within  a  segment  -  a  <  :  <  a  and  zero  elsew’here  then 
B^(k)  =  sin  n^a/n^a  and  so 


In  general  the  sum  in  (148)  can  be  written  as  an  integral,  because 


B  (k) 


J  p^(c)  exP  ~  df 


and  so 


I  <->"V> 

nk=-“ 


exp  iVJ\<V 


O 


j  p^(f')dv  (-)  k  exp  in^(6  -  (JlJn^Cr^C) 

V"“ 

77 

j  Pk(f)  exp[-  ir^C  sin(6  -  d  0 


(154) 


and  this  leads  to  the  third  and  last  example  in  this  section,  the  case  of  a  Gaussian 
step  angle  density  with  small  variance.  Let  P  (l?)  be  the  theta  density  of  section  4  (99). 
Then 

j  P^f;)  exp  -  ir^'  s  in ( •■  -  :)d: 


t 


exp 


(c  -  2nr)  ‘ 


exp[-  ir^S  sin  (e  - 


f  exn  -  -~p  exp[-  ir  C  sin(B  -  l>)]dc  .  (155) 

■  2o“  k 


It  will  now  he  assumed  that  the  deiisi t  v  function  is  sharply  peaked  around  $  *  0  ,  ie  c 
is  small  so  that  pi  1  will  also  he  sharply  peaked  around  e=0  and  will  thus  be  very 
small  outside  the  positive  halt  plane.  Hence  making  the  approximation 
sin("  -  :  1  ■»  siti  -  ■'  c»»s 


30 


/ 


P,  (4)  exp  -  ir,  C  sin  ('•  -  4)d;  ~  — = 

k  r  k  o/H 


<P(~irk 'sinh)j  exp-^-j  exp(irkii 

—CO 

!4i  • 


cos  :)d: 


exp  -  ir  sin  t?  exp<-  r  a  cos 


{- 


2  2 


0  56) 


And  so 


p(6) 


I  £i  r-7  |exp(-  isr  sin  v)  exp  -  t 2c2  cos2e  ~  |  , 


057) 


where  s  = 


N  ,  N 

y  r  and  t  =  V  r 
T  k  •  u 


The  integration  is  straightforward  and  leads  to: 


,  - ,  Is  I 

pO)  -  —  — — rye 

*  -  cos  6 


xp/-  tan“9  )  , 

\  2t~a‘ ' 


0  58) 


which  for  6  small  gives: 


p(-) 


-k £  • 


059) 


Note  that  in  (158)  there  is  an  approximately  Gaussian  function  centred  on  5=0  and  a 
second  centred  on  6  =  ±ir  .  The  second  is  positive  whereas  according  to  the  analysis  of 
this  section  it  should  be  negative.  The  reason  for  this  is  that  in  making  the  approxima¬ 
tion  sin  $  21  $  and  cos  $  —  1  the  periodic  properties  of  sin  4  and  cos  4  are 
destroyed.  Near  i  =  ±ir  ,  sin  $  =*  -4  and  cos  4  =  -1  and  if  this  is  inserted  into  (156) 
it  will  be  found  that  -p(9)  results  in  the  negative  half  plane. 


5 . 3  The  phase  density  function  on  the  whole  complex  plane 


In  order  to  obtain  p(6)  over  the  whole  plane  it  is  necessary  to  integrate  (130) 
and  this  seems  to  be,  in  general,  intractable.  One  method  is  to  note  that  the  random 
walk  has  a  finite  span,  and  p(a,0)  =  0  for  a>jr|  +...+r^|.  Then 


4  [l  I  rk  cos 
-  I  rk  t0 

u,  =  0  or 


\\  +  l  rkcos  can  be  expanded  as  a  Fourier  series  on  the  interval 

+  1  r^  ,  and  the  result  will  be  valid  for  all  u^  with  the  sole  exception 
wr  where  the  Fourier  series  will  give  J  J  r^  instead  of  I  rk  • 


of 


Alternatively,  an  integral  representation  could  be  used.  The  technique  of  expanding  a 
function  on  the  span  of  a  random  walk  has  been  employed  before  bv,  for  example.  Slater*2. 
The  resulting  series  for  p(9)  will  not  be  given  here;  it  seems  to  be  beyond  all  hope 
of  summation  and  practical  usefulness.  As  pointed  out  in  the  Introduction,  the  main 
interest  is  in  computing  the  first  order  term  in  p(6)  for  a  small  deviation  from 
uniformity  in  the  step  density.  As  shown  in  section  4,  for  a  Gaussian  step  angle 
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density  (the  practical  case  of  interest)  the  step  angle  density  can  be  very  accurately 

described  by  the  constant  and  first  harmonic  (fundamental)  in  the  Fourier  expansion  of 

2 

P  ( ; ) ,  provided  that  o“  >  2r  .  Fortunately,  to  obtain  p(0)  in  this  case  it  is  not 
necessary  to  integrate  (130),  and  it  will  be  shown  that  the  behaviour  of  p(t’)  can  be 
deduced  from  the  odd-numbered  harmonics  alone,  from  equation  (148).  Substituting 
(1  +  b  cos  f ) / 2 n  for  P($)  (and  thus  assuming  equal  length  steps)  into  equation  (148) 
leads  to  an  awkward  analysis,  however,  and  a  more  elegant  method  for  obtaining  p(i-)  for 
small  non-uniformity  is  the  following: 

Assume  that  each  step  has  equal  (and  arbitrary)  length.  This  simplification  is  not 
strictly  necessary  but  gives  tidier  expressions,  and  is  easily  removed.  Let  the  angle 
probability  density  function*  of  each  step  be 

P^ff)  =  exp(a  cos  ip)  /2vT^(a)  ( 160) 


where  I^fa)  is  t*ie  modified  Bessel  function  of  the  first  kind  (order  zero).  Then 

IT 

P  (<)>)  >0  and  /  P  (4)  =  I  . 

-IT 

Note  that  the  probability  density  functions  for  each  step  are  'in  phase'.  The  corres¬ 
ponding  Fourier  series  is 


V*> 


■  i  £»„ 


exp  inf 


where  B 


li 

2*In(a)  I 


exp(a  cos  f  -  inf)df 


and  so 

Bn  =  In(a)/I0(a) 

Hence 


Pk(4) 


1  (a) 

yir exp  in* 


(161) 


(162) 


So  in  equation  (148)  it  is  necessary  to  find  the  sum 


El  (a) 

<->"  yiy 


S  = 


exp  inP 


This  is  easily  found  by  Neumann's  addition  theorem  (Watson**’*  §11.2  (1)). 

Since  m 

C"» 

(z)  exp  inf 


JnM  = 


ZJ  (Z)J  ( 
n  n 


(163) 


(164) 


*  This  function  has  been  plotted  in  Fig  8  for  a  =  0  to  1.0  in  steps  of  0.1.  This 
function  is  Von  Mises  density  function  (Ref  4,  equation  (12), also  Ref  19,  section  3.4.9, 
page  57) . 


n — 

uJ  =  ^  Z  ♦ 


where  w  =  JZ~  *  z~  -  2Zz  cos  <f  and  Z,z  have  general  complex  values,  putting  7.  =  ia  , 
z  «  5  ,  $  «  0  + n/2  gives 


(r^  2 

Jr,~  -  a  +  2ia.r,  sin  e 


( 185) 


Hence  if  q ( 0 )  is  the  'odd  numbered  harmonic  content'  of  p(B)  ,  then 

,N 


r  .  t 

LVa)J 

1  7.  30  [ 

a~  +  2ia£  sin 


.  (188) 

It  is  now  assumed  that  (i)  a  is  small  <0.1  say  and  (ii)  N  is  large  >10  say. 

Under  these  conditions  it  can  be  shown  (see  Watson  §13.48,  page  421)  that  the  signifi- 

n  2 

cant  part  of  the  integrand  is  that  for  which  -  a  +  ?iar  sin  is  small,  and  then 

N-l 


£j0(Jr,2  -  a2  *  2ia-  sin  ')j  J ,  ( J ' 2  -  a2  ♦  2ia-  sin  o)^2"- 

(N  -  J) 


a  +  :ia.1  sin 


!  exP 


[-* 


(r.“  -  a"  +  2iaf  sin 


( 187) 


In  connection  with  this  result  it  is  easy  to  see  that  for  ' z\  small  J^(z)  "  1  “  ~ 


hence  [V2^  '  ~  t)  - 


So 


2q  (0) 


if  I  ~|  Na  cos  0  f 

27[yiiJ  —2—1 


exp 


-  1)  ^  2 

■>  ;  (C  -  a  +  2  i  a  f.  rfin  ft)dt  .  (168) 


Also,  for  small  a 


r0(a)  -  1  +  a  /4  . 


(189) 


Hence 


£lQ(a)J  ~  (1  +  a“/4)N  ~  exp(NaV4)  a  small 


(170) 


And  so,  if  J  is  neglected  in  comparison  with  N  : 


.  /a\  1  Na  cos  6  f  N  ,r2  .  .  .... 

2q(0)  —  - - -  I  exp  -  j  (r.  +  2iaf  sin  'Md’.  , 


(171) 
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a  cos  8 
4 


f N  /  Na2  sin2e  \ 

J-exV(-  - - - j  , 


and  this  is  an  interesting  result. 

So  far  it  has  been  assumed  that  N  is  large  and  a  is  small.  Now  suppose  that  Na' 
is  large,  and  then: 


t/i  ‘4  »£)  • 


Note  that  this  is  a  '  convergent  sequence,  hence 


CXj 

-i£ 


(”)■-('?  -  nr)  , 


which  gives  delta  functions  of  amplitude  +1  at  •'<  =  0  and  -1  at  ?  s  -n 

Since  p(")  >  0  the  even  harmonics  must  supply  a  delta  of  amplitude  at  least  +1  at 

v  -  -t  ,  and  hence  also  a  delta  of  amplitude  at  least  +1  at  0  =  0.  Since 
’! 

I  p(  )d0  =  I  the  amplitudes  of  the  deltas  at  0=0  and  -  are  each  +1  ,  hence 


as  N  -»  * 


£ 1 


!(■■  -  2n+) 


Therefore,  for  a  small  departure  from  uniformity  no  matter  how  small  the  phase  probability 
density  function  of  the  resultant  tends  towards  a  delta  function  as  the  number  of  steps 
increases  without  limit.  When  the  phase  density  function  of  each  step  are  ’in  phase’  the 
effect  of  the  random  walk  is  to  amplify  the  non-uniformity  in  the  density  of  each  step. 

The  first  harmonic  in  qO)  will  next  be  obtained. 


a  (n  „ 

-  zj~  r- 


Na'  sin'' 
exp - t - - 


=  £  (i  V  L  ad)  jl 

L  \  4/k! 


cos  v  sin 


and  also 


L_(i-  lit-  it  ,  _Jl_  cos  ,p  +  \ 

: or (k  +1)12  k  +  i  s  '  "7 


see  Whittaker  and  Watson 


,  page  191,  example  13. 
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Hence 


.  2k,  2r(k  +  })  <  1  k 

sin  9  *  rlj')F(k~"+  0  { I  ~  lT~i~T  cos  29 


re 


,  ■  2k 

cos  8  sin  8 


Hence 


rllW'MT  {COS  9  "  iTTT  cos  e  cos  28  +  . . .  j 
r([)r(k  +\)  CQS  9  +  hi8her  harmonics 


oo 

q(e)  -  f/?cos  9  Y,  ("T")  ^T7(1T^"2T  +  higher  hannonics 

k=0 

=  ^ J ~  cos  0  J»2,  -  j  +  higher  harmonics.  » 


where  F.  is  a  confluent  hypergeometric  function  (Whittaker  and  Watson  ,  §16.1, 

2 

Evidently,  for  small  Na  ,  we  have 

q(0)  £jl  cos  , 

2 

p(0)  =  ~  £l  +  cos  e{*  ~  +  •••}  +  higher  harmonics  j  . 

2 

Hence  the  behaviour  of  p(8)  is  given  by  (184)  for  small  values  of  Na  (say  for 


Na  <  1 )  and  by 

,„N  a  f¥  Na^e^ 

P(0)  =  jJ-  exp  -  —3—  , 

2  2 

for  large  values  of  Na  (Na  >  10).  In  each  case  N  is  large  and  a  small. 
The  step  phase  density  is 


i  vIn<a) 

=  27  LT^Texp  ln*  • 


V*> 


i  r  2va)  i 

“  27['  +  T^TaT  cos  ♦  +  •••]  * 


since  I  (a)  =  T  (a).  For  small  a 
-n  n 


and  so 


V*>  „  . 
2 


Pk(<f)  ~  fl  +  a  cos  4>]  +  0(a; 


(179) 

(180) 

(181) 

(182) 

page  338) . 

(183) 

(184) 

(185) 

(162) 

(186) 

(187) 

(188) 
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ror  the  theta  density  of  section  4  for  a  >  2  it  : 


Pk(4)  "  j^l  +  2  exp  -  —  cos  <fj  +  0(exp  -  a")  .  (184) 

And  so  the  results  above  are  valid  for  a  theta  density  with  large  variance  if  we  put 

2 

a  =  2  exp  -  ~  .  (140) 

The  analysis  presented  here  is  not  particularly  rigorous  and  so  it  is  checked  in  the 
next  section  by  comparing  results  (184)  and  (185)  against  the  phase  probability  density 
function  obtained  by  numerically  computing  random  walks  with  a  Gaussian  step  phase 
density. 

6  NUMERICAL  TESTS 

In  principle  it  is  a  straightforward  matter  to  numerically  compute  the  phase 
probability  density  of  the  resultant  of  a  non-i sotropic  random  walk  using  a  Monte-Carlo 
technique.  Suppose  that  it  is  desired  to  compute  the  phase  density  of  an  N-step  walk 
with  theta  phase  density.  A  quasi-random  number  generator  with  a  Gaussian  density  sup¬ 
plies  the  phase  angle  of  each  step.  The  fact  that  one  computes  the  changes  in  x  and  v 
coordinates  of  the  resultant  (with  the  addition  of  each  sten)  from  Ax  =  r  cos  f-  and 
Ay  =  r  sin  0  ensures  that  the  tails  of  the  Gaussian  are  wrapped  back  on  to  the  interval 
because  of  the  periodicity  of  sin  0  and  cos  ?  .  The  theta  probability  density 
is  therefore  generated.  When  N  steps  have  been  added  together,  the  phase  of  the 
resuLtant  is  computed  and  stored.  This  process  is  repeated  many  times  and  when  suf¬ 
ficient  walks  have  been  computed  a  histogram  .s  constructed  from  the  stored  phases.  This 
histogram  is  then  an  approximation  over  finite  intervals  of  the  probability  density  of 
the  phase  of  the  resultant  of  the  walk.  This  is  the  'Monte-Carlo'  method. 

In  this  section  histograms  generated  using  the  above  method  are  compared  with  the 
analytic  results  of  the  previous  section.  There  are  certain  numerical  difficulties 
which  should  be  mentioned.  In  particular  the  quasi-random  number  generators  found  in 
the  libraries  of  most  computers,  and  used  to  generate  the  numbers,  are  usually  not 
entirely  satisfactory.  This  is  despite  the  fact  that  these  generators  pass  the  standard 
tests  for  'randomness'  of  the  numbers  generated.  In  practice  the  higher  order  moments 
of  these  numbers  are  far  from  ideal  and  if  a  large  number  of  random  numbers  is  required 
then  correlations  between  the  numbers  distort  the  random  walk  results.  Great  care  is 
necessary.  It  was  found  necessary  to  re-randomise  the  output  of  the  particular  generator 
used  by  accepting  or  rejecting  the  numbers  according  to  an  independent  table  of  Os  and  Is. 
This  is  easily  constructed  by  tossing  a  coin.  The  length  of  such  a  table  is  important. 

It  should  not  be  too  long  nor  too  short  and  experiments  must  be  carried  out  to  determine 
the  optimum.  In  this  case  about  100  binary  values  proved  optimum.  Greater  lengths 
produced  greater  smoothing  of  the  probability  density  of  the  numbers,  but  introduced 
skewness.  The  random  number  generator  used  produced  numbers  uniformly  distributed  on 
5  the  interval  [0,lj.  Gaussian  variables  were  produced  from  this  by  the  method  of  Box 

J  g 

and  Muller  ,  the  so-called  'direct  method'. 
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In  the  computations  which  follow  the  histograms  were  each  constructed  from  4096 
walks.  Each  walk  consists  of  from  1  to  10000  steps.  Thus  in  the  most  extreme  case  about 
4  x  ,07  random  numbers  were  required.  Each  histogram  has  either  32  or  64  sub-intervals 
covering  the  interval  [-utnj.  Thus,  on  the  average,  each  sub-interval  contains  the  phase 
of  the  resultant  of  128  or  64  walks.  If  the  number  of  sub-intervals  is  equal  to  some 
power  of  2  this  makes  the  computation  of  the  harmonic  content  easy  via  a  discrete  (fast) 
Fourier  transform. 


The  objective  of  the  numerical  tests  is  to  check  the  validity  of  formulae  (159), 
(184)  and  (185)  and  they  are  treated  in  that  order. 

6. 1  Gaussian  step  density  of  small  variance 

For  equal  length  steps  formula  (159)  gives: 


p(6) 


(191) 


where  a  is  the  variance  of  the  phase  distribution  of  each  step,  and  is  assumed  to  be 

small  enough  so  that  most  of  the  Gaussian  distribution  is  contained  within  *^/2  .  A 
2  2 

value  of  (n/6)'  was  chosen  for  o  and  the  resulting  phase  density  of  the  resultant, 
was  computed  for  I,  2,  5,  10,  20  and  50  steps.  The  case  of  1  step  was  chosen  as  a  check 
on  the  random  walk  generator.  It  will  be  noted  that  (159)  is  valid  for  any  number  of 
steps.  The  results  are  shown  in  Figs  9  to  14.  In  each  case  the  histogram  is  computed  hv 
the  Monte-Carlo  method  and  the  super-imposed  Gaussian  curve  is  that  given  by  (1Q1)  with 
the  appropriate  value  of  N  .  It  will  be  seen  that  in  each  case  formula  (]9l)  fits  the 
computed  histograms  quite  closely.  Note  that,  for  a  constant  value  of  the  proba¬ 

bility  density  becomes  'sharper'  as  N  increases. 


Table  I  lists  the  variance  of  the  phase  density  of  the  resultant  comnuted  from  the 
histogram  and  also  that  given  by  (191).  Again  there  is  close  agreement.  The  proportional 
error  becomes  larger  as  N  increases  and  is  caused  by  the  finite  and  fixed  number  of 
sub-intervals  in  the  histogram. 


Table  I 


Number  of 
steps 

1 

o'  from  histogram 

o'  f rom  (  1  8  | ) 

1 

0.270 

0.274 

2 

0.136 

0.  1  37 

5 

0.0542 

0.0548 

10 

0.0286 

0.0274 

20 

0.0145 

0.01 37 

50 

0.00645 

0.00548 

In  all  cases  the  step  phase  density  has  a  variance  of  (r/63 
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6 . 2  Gaussian  step  density  of  large  variance  (Na"  <  I) 


When  the  underlying  step  phase  density  has  a  large  variance  and  the  number  of  steps 

2  2  2 

N  is  such  that  Na  <  I  where  a  =  2  exp(-a  /2)  ,  a  being  the  step  phase  variance, 
then  the  phase  density  of  the  resultant  is  given  by  (184) 
ie 


p(Q) 


i  r 

,  a  / 

Na 2  \  , .  ,  .  1 

L 

1  +  y  vNtt  cos  P  1 

1 - j-g-  +  . .  •  J  +  higher  harmonics  1 

(184) 


This  formula  has  been  tested  by  generating  random  walks  having  2,  4,  6,  8,  10,  15  and  20 

2 

steps  for  values  of  Na  of  1,  0.1,  and  0.01;  a  total  of  21  cases.  Although  the  analysis 

is  only  valid  for  large  N  the  computed  histograms  show  remarkable  agreement  with  (184) 
for  only  2  steps.  It  was  therefore  decided  to  include  walks  with  a  small  number  of  steps 
as  well  as  a  moderate  number:  it  might  be  expected  that  N  =  10  would  be  the  lower  limit 
on  the  validity  of  (184). 

2 

Figs  15,  16  and  17  show  the  histograms  obtained  for  N  =  2,  10  and  20  and  Na  =1. 

2 

Figs  18,  19  and  20  show  histograms  for  N  =  2,  10  and  20  for  Na  =0.1  and  Figs  21,  22 

2 

and  23  for  N  =  2,  10  and  20,  Na  =  0.01 .  In  each  case  the  superimposed  constant  plus 
first  harmonic  is  that  obtained  from  the  histogram  by  a  discrete  Fourier  transform.  It 
represents  the  'filtered'  histogram  with  all  frequencies  greater  than  the  fundamental 
cut  off. 


Table  2  lists  the  amplitude  of  the  first  harmonic  for  each  computed  case.  Those 
cases  marked  with  a  star  are  those  in  Figs  15  to  24. 


Table  2 


Number  of 
steps 

Na2 

1.0 

0.  1 

0.01 

2* 

0.904 

0.301 

0.109 

4 

0.861 

0.277 

0.080 

6 

0.851 

0.288 

0.105 

8 

0.844 

0.282 

0.100 

10* 

0.844 

0.300 

0.099 

15 

0.840 

0.282 

0.101 

20* 

0.832 

0.278 

0.084 

CO 


Values  of  the  fundamental  (first  harmonic)  given  by  (184)  are  0.831,  0.278  and 
2 

0.0886  for  Na  =  1  ,  0.1  and  0.01  respectively.  Note  that  the  amplitude  of  the  first 

a  r—  /  Na  2  \ 

harmonic  is  taken  to  be  —  I  I  -  — +  •••  1  and  the  factor  of  Jit  has  not  been 
included . 
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There  is  reasonable  agreement  between  (184)  and  the  Monte-Carlo  computations,  even 

2 

for  small  values  of  N  .  It  will  be  noted  that  at  Na  =1  there  is  also  a  significant 

second  harmonic  component,  of  amplitude  approximately  1/5  of  the  first  harmonic. 

2 

Na  =  1  ,  therefore,  represents  the  approximate  upper  limit  of  the  validity  of  (184). 

2 

6 . 3  Gaussian  step  density  of  large  variance  (Na  >  10) 

In  this  case  we  have  formula  (185) 

a  fN  Na2S2 

p(0)  ~  jJ-  exp - —  ,  (185) 

and  this  represents  the  second  of  the  two  asymptotic  limits  of  p (6 )  when  the  under¬ 
lying  step  phase  density  has  a  large  variance.  It  has  been  tested  by  computing  random 
walks  for  Na2  =10,  N  =  10,  102  and  103  (Figs  24,  25  and  26),  for  Na2  =  20,  N  =  102, 

103  and  104  (Figs  27,  28  and  29),  and  for  Na2  =  100,  N  =  102,  103  and  104  (Figs  30,  31 

and  32).  In  each  case  the  superimposed  Gaussian  density  function  is  obtained  from  ()85) 

2 

with  the  appropriate  value  of  Na  .  There  is  reasonable  agreement  in  all  cases.  Table  3 

lists  the  variance  obtained  from  the  histograms.  The  variance  according  to  (185)  should 

2 

be  0.2,  0.1  and  0.02  for  Na~  =  10,  20  and  100  respectively.  There  is  thus  reasonable 

•  •  •  2 
quantitative  agreement  which  seems  to  be  better  for  the  larger  values  of  Na",  as  is 

4 

to  be  expected.  The  values  for  N  =  10  are  probably  distorted  by  imperfections  in  the 
random  number  generator  which  show  up  for  large  numbers  of  steps. 

Table  3 


Number  of 
steps 

Na2 

10 

20 

100 

10 

0.207 

- 

- 

102 

0.273 

0.112 

0.020 

103 

0.270 

0. 1  10 

0.021 

I04 

. 

0.084 

0.018 

7  CONCLUSIONS 

It  has  been  demonstrated  that  a  theory  of  the  phase  probability  densitv  function  of 
the  resultant  of  a  non-isotropic  two-dimensional  random  walk  is  possible  without  recourse 
to  the  central  limit  theorem.  Although  the  most  general  case  seems  to  be  totally  intract¬ 
able,  in  the  case  of  a  Gaussian  step  phase  probability  density  function  (theta  densitv) 
the  following  results  have  been  obtained: 

(I)  If  the  variance  of  the  phase  probability  density  of  each  step  (a2)  is  small, 

1  <  n/3  say,  then  for  N  equal  steps  the  phase  probability  densitv  of  the 
resultant  is  giver  by: 
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i  .  'n  /  N'-r’\ 

rn  o  exp(“  2o2j 


l 


and  this  is  valid  for  any  number  of  steps. 

2  9 

(2)  (a)  If  the  variance  o  of  the  step  phase  density  is  large  o~  >  2 t  say,  and  the 

2  2 

number  of  steps  N  is  such  that  Na"  <  1  where  a  =  2  exp(-o“/2)  then  the  phase 
probability  density  of  the  resultant  is  given  by 


27  ['  +  f  ^  cos  e  (!  "  TT  +  •••)] 


This  formula,  derived  for  N  large,  is  also  reasonably  good  for  N  small  (*2). 

? 

(b)  If  N  is  such  that  Na“  >  10  then 


a  fN  /  Na292  \ 

p(,)  - 


and  this  is  valid  for  large  N  . 

These  three  cases  are  of  physical  interest  and  are  useful  results. 

In  conclusion  it  is  pointed  out  that  if  Pearson's  drunk"  is  slightly  biassed  in  his 
direction  of  stagger  (perhaps  he  is  trying  to  make  for  home)  then  if  he  takes  enough  steps 
it  is  very  probable  that  he  will  go  in  the  desired  direction,  provided  he  does  not  collapse 


of  course. 
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